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Abstract 

A key aspect of fluid dynamics is the correct definition of the phase-space Lagrangian dynamics 
which characterizes arbitrary fluid elements of an incompressible fluid. Apart being an unsolved 
theoretical problem of fundamental importance, the issue is relevant to exhibit the connection 
between fluid dynamics and the classical dynamical systems underlying incompressible and non- 
isothermal fluid, typically founded either on: a) a configuration- space Lagrangian description of 
the dynamics of fluid elements; b) a kinetic description of the molecular dynamics, based on a 
discrete representation of the fluid. The goal of this paper is to show that the exact Lagrangian 
dynamics can be established based on the inverse kinetic theory (IKT) for incompressible fluids 
recently pointed out (Ellero et al., 2004-2006, Q]). The result is reached by adopting an IKT 
approach based on a restricted phase-space representation of the fluid, in which the configuration 
space coincides with the physical fluid domain. The result appears of potential importance in 
applied fluid dynamics and CFD. 

PACS numbers: 47.10.ad,05.20.Dd 
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I. INTRODUCTION 



As it is well known, the description of fluids can be performed choosing either an Eulerian 
or a Lagrangian point of view. The two approaches are equivalent and, if fluid dynamics 
were fully understood, one should be able to translate Eulerian properties into Lagrangian 
ones and viceversa. At least for the treatment of turbulent fluids, we are still quite far from 
this point The main historical reason of this situation can be understood by looking at 
the customary statistical approach based on so-called velocity probability density function 
(pdf)- method for an incompressible fluid (for a review see for example [^j]). In fact - as 
it is well known - in this approach the time evolution of the pdf which advances in time 
the average (in a stochastic sense) fluid velocity is determined by a Fokker-Planck transport 
equation. As a consequence the corresponding Lagrangian characteristics are necessarily 
stochastic in nature and therefore difficult to handle. This explains why in the literature, 
for this purpose, stochastic models of various nature have been adopted, which, however, 
typically rely on experimental verification rather than (uniquely) on first principles. 

Recently, however, an important breakthrough has been achieved by the discovery of the 
so-called inverse kinetic theory (IKT) approach for incompressible fluids 3), \a, [f], Q] which 
permits a straightforward connection between Eulerian and Lagrangian descriptions. This 
is achieved by identifying the relevant fluid fields, which are assumed to be defined in a 
suitable domain Q C IR 3 (fluid domain), with appropriate moments of a suitably-defined 
kinetic distribution function (KDF) /(x,t) [with x = (r, v) e T, x and T denoting a suitable 
state-vector and an appropriate phase-space] which is assumed to advance in time by means 
of Vlasov-type kinetic equation. In such a case, the time-evolution of the KDF is determined 
by a kinetic equation which, written in the Eulerian form, reads 

L/(x,t) = 0. (1) 

Here /(x, t) denotes the Eulerian representation of the KDF, L is the streaming operator 
Lf = J^/ + ^ ■ {X(x, t)f}, X(x, t) = {v, F(x, t)} a suitably smooth vector field, while 
v and F(x, t) denote respectively appropriate velocity and acceleration fields. As a main 
consequence the approach can in principle be used to determine in a rigorous way the 
Lagrangian formulation for arbitrary complex fluids. Although the choice of the phase 
space T is in principle arbitrary, in the case of incompressible isothermal fluids, it is found 
4| that the phase-space T can always be reduced to the direct-product space T = Q x V 



(restricted phase-space), where Q, V C M. 3 , Q is an open set denoted as configuration space 
of the fluid (fluid domain) and V is the velocity space. This type of approach (based on a 
restricted phase-space IKT formulation) will be adopted also in the sequel. 

The main motivation [of this work] is that some of the general understanding recently 
achieved in simple flows by means of the IKT approach could also give a significant con- 
tribution to a wider range of problems. In the sequel, we will concentrate on the issue of 
a consistent formulation for fluid dynamics based on a phase-space (IKT) description of 
incompressible fluids, whereby its pressure, velocity (and possibly also thermal) fluctuations 
are consistently taken into account. In particular we intend to show that the theory ad- 
mits a well-defined Lagrangian formulation. In particular, the corresponding Lagrangian 
trajectories can be interpreted in terms of the phase-space dynamics of suitable, classical 
molecules, i.e., point particles whose dynamics is determined by the phase-space Lagrangian 
characteristics. The motion of these particles (rather than that of the fluid elements), as 
they are pushed along erratic trajectories by fluctuations of the fluid-field gradients (in par- 
ticular, characterizing the fluid pressure and temperature), is fundamental to transport and 
mixing processes in fluids. In this regard, it is well known that the interaction between 
(deterministic and/or turbulent) fluctuations of the fluid fields and fluid particles still es- 
capes a consistent theoretical description. Being a subject of major importance for many 
environmental, geophysical and industrial applications, the issue deserves a careful investiga- 
tion. A key aspect of fluid dynamics is, therefore, the correct definition of the (phase-space) 
Lagrangian dynamics which characterizes incompressible fluids. The customary approach 
to the Lagrangian formulation is based typically on a configuration-space description, i.e., 
on the introduction of the configuration-space Lagrangian characteristics r(t), spanning the 
fluid domain Q. Here r(t) denotes the solution of the initial-value problem: 



with r Q an arbitrary vector belonging to the closure Q of Q and V(r, t) being the velocity 
fluid field, to be assumed continuous in Q and suitably smooth in Q. The purpose of this 
paper is to achieve a phase-space Lagrangian formulation for incompressible thermofluids, 
extending the formulation previously developed, adopting, instead, an extended-phase-space 
formulation [8j . The present formulation permits to advance in time the relevant fluid fields 




(2) 
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by means of the set of Lagrangian equations defined by the vector field X(x, t), namely 



where x Q is an arbitrary initial state of T (closure of the phase-space T). The goal of this 
paper is to show, in particular, that the exact (phase-space) Lagrangian dynamics can be 
established based on the inverse kinetic theory (IKT) for incompressible fluids in such a way 
that the phase space T coincides with the direct product space T = Q x V, Q being the fluid 
domain and V (velocity space) the set R 3 . The result appears relevant in particular for the 
following reasons: 1) the Lagrangian dynamics here determined permits to advance in time 
self-consistently the fluid fields, i.e., in such a way that they satisfy identically the requires 
set of fluid equations. For isothermal fluids, this conclusion is consistent with the results 
indicated elsewhere by Tessarotto et al. [5|; 2) the Lagrangian dynamics takes into account 
the specific form of the phase-space distribution function which advances in time the fluid 
fields; 3) the theory permits an exact description of the motion of those particles immersed 
in the fluid which follow the Lagrangian dynamics (classical molecules). In detail the plan 
of the paper is as follows. First, in Section 2, an IKT for incompressible thermofluids, 
adopting a restricted phase-space description (analogous to that developed in Ref. [5]), is 
presented. This permits us to determine the appropriate form of the vector field X(x, t). For 
the sake of illustration, we shall consider in particular the case in which the KDF is identified 
with a local Maxwellian distribution function. Second, in Sec. 3, the new set of Lagrangian 
equations is presented, which are proven to advance uniquely in time the relevant fluid fields 
of an incompressible thermofluid. 

II. EULERIAN FORMULATION OF IKT 

A first key issue is to prove that the IKT approach developed for isothermal incompressible 
fluids can also be achieved by adopting the same type of restricted phase- 

space formulation, obtained identifying the phase space with T = Q x V, where Q coincides 
with the fluid domain and V = M 3 is a 3D velocity space. For definiteness, here we shall 
assume that the relevant fluids for a thermofluid, namely {p = p Q > 0, V,p > 0, T > 0, St} , 
i.e., respectively the constant mass density, the fluid velocity, pressure, temperature and 




(3) 
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entropy, in the open set satisfy the so-called non-isentropic and incompressible Navier- 
Stokes- Fourier equations (INSFE), i.e., 



V- V = 0, 




Po 



1 



[Vp - f] + z/V 2 V, 



(4) 
(5) 




) 



2 



1 



2c p \dx k dxi 



+ 



J = K, 



(6) 



PoC; 



V 



(7) 



where j^V is the fluid acceleration, 757 = ^ + V-V the convective derivative and St a 
functional, to be suitably defined, of an appropriate set of fluid fields. These equations 
are assumed to satisfy a suitable initial-boundary value problem (INSFE problem) so that 
a smooth (strong) solution exists for the fluid fields {p = p Q > o, V,p > 0, T > 0} . Here 
the notation is standard. Thus, Eqs. (HJ), (JSj) and (EJ) denote respectively the so-called iso- 
choricity condition, the forced Navier-Stokes equation, written in the Boussinesq approxi- 
mation, and the Fourier equation. As a consequence, in such a case the force density f reads 
f = p g (1 — k p T) + Fi, where the first term represents the (temperature-dependent) grav- 
itational force density, while the second one (fi) the action of a possible non-gravitational 
externally-produced force. Hence f can be written also as f = — V0 + {r, where <fi = p a gz 
and f# = — p gA;pT + fi denote respectively the gravitational potential (hydrostatic pressure) 
and the non-potential force density. Moreover, in Eq. ([6]) J is the quantity of heat generated 
by external sources per unit volume and unit time. Finally, Eq.([7]) defines the so-called 2nd 
principle for the thermodynamic entropy St- For its validity in the sequel we shall assume 
that there results everywhere in Q x / 



(externally heated thermofluid) . In these equations g,k p , u,x and c p are all real constants 
which denote respectively the local acceleration of gravity, the density thermal-dilatation 
coefficient, the kinematic viscosity, the thermometric conductivity and the specific heat at 
constant pressure. Thus, by taking the divergence of the N-S equation ([5]), there it follows 
the Poisson equation for the fluid pressure p which reads 




(8) 



V 2 p = -p V ■ (V • VV) + V ■ f , 



(9) 
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with p to be assumed non negative and bounded in Q x I. In this Section we shall assume 
that /(x,£) is a solution of the Eulerian kinetic equation ([I]) defined in a suitable extended 
phase-space T x I, where / C R is a suitable time interval. In such a case, we intend to 
show that /(x,i) can be defined in such a way that the fluid fields V,pi can be identified 
with its velocity moments J dv(j(x,t)/(x,£), where respectively G(x,t) = l,v, u 2 /3 and pi 
is the kinetic pressure to be defined as: 

Pi =Po(t) + p- (p + poT/m. (10) 



Here po(t) (to be denotes as pseudo-pressure {§]) is an arbitrary strictly positive and suitably 

smooth function defined in I. Moreover, m > is a constant mass, whose value remains in 

principle arbitrary. In particular it can be identified with the average mass of the molecules 

forming the fluid. Finally, the thermodynamic entropy St can be identified with the Shannon 

statistical entropy functional S(f(x,t)) = — f dx/(x,£) ln/(x,t), provided the function po(t) 

r 

is a suitably prescribed function and /(x,i) is strictly positive in the whole set T x /. To 
reach the proof, let us first show that, by suitable definition of the vector field F(x, t), a 
particular solution of the IKE (pQ) is delivered by the Maxwellian distribution function: 

/ M (x,t)= P ° 3 expJ — ^— I. (11) 

Here u = v — V(r, t) and Vth, P = ^/2j?i(r, t)/p are respectively the relative and the thermal 
velocities. The following theorem can immediately be proven: 

Theorem - Restricted phase-space INSFE-IKT Let us assume that: 1 ) the INSFE 
problem admits a smooth strong solution in Y x I, such that the inequality (TJj) is fulfilled; 2) 
the vector field F is defined as 

F(x,f;/)=F + Fi, (12) 

where F , F 1 read respectively 

F (x,t; /) = - [V-S - V Pl + f fl ] + u • VV+z/V 2 V, (13) 

Po 

Px(x,t; /) = iu {-A+lv • Q - i [V-HJ -q) + ^V-H - |j , (14) 
2 pi Pi — J 2pi — K 2 J 



D 1 



d 

A = -^(po+p)-y- 



— V-— f R -zA7 2 V 



+ = TrrPi- (15) 

m Dt 
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3) in O x / the KDF /(x,£) admits the velocity moments G(x,t) = l,v,-u 2 /3,u^ and 
uu; i/ius, we denote in particular Q = J d 3 vu^f and II = f ££ 3 i>uu/; ^ £/ie iCDF /(x,£) 
satisfies identically in in Q x I the constraint equation: 

p= I dv/(x,t) = p > 0, (16) 



where p is a positive constant; 5) the entropy integral S(f(x.,to)) = — J dx/(x,t ) ln/(x,i ) 

r 

exists; 6) in the time interval I the pseudo-pressure Po(t) is defined so that there results 
identically: 

/ dr- _ Pl + V-Q— VpQ =0. (17) 
Jn Pi l dt P 

It follows that : A ) the local Maxwellian distribution function 177]) is a particular solution 
of the IKE (Qp if an only if the fluid fields {p = p D > 0, V,p, T} satisfy the fluid equations 
0)-(0). In snc/i a case there results identically Q = 0, II = 0; B) for an arbitrary KDF 
/(x, t), which satisfies identically assumptions 4)-6), the velocity-moment equations obtained 
by taking the weighted velocity integrals of Eq. (QP with the weights G(x, t) = 1, v, u 2 /3 deliver 
identically the same fluid equations ([^-(ED/ C) the Shannon entropy S (/) is a monotonic 
function of time, i.e., there results (\/t El)'- 

p(f)>° ( 18 ) 



(H-theorem) , where there results in particular for isothermal fluids (see also [8]) J^S (/) = 
(constant H-theorem). This permits us to identify St = S (/) so that the 2nd principle 
[i.e., the inequality 0/ is satisfied too. 

PROOF First, let us assume that a strong solution of the INSFE problem exists which 
satisfies identically Eqs.(jl])-([6]) in the set Q x I. In such a case it is immediate to prove 
that /w(x, t) is a particular solution of the inverse kinetic equation (Tj[|). The proof follows 
upon invoking Eqs. (|T2|) - (fl5|) for the vector field F(x, t; f) and by direct substitution of the 
distribution /iw(x, t) in the same equation (Proposition A). Instead, if we assume that in 
r x I, f = /a/(x, t) is a particular solution of the inverse kinetic equation, which fulfills 
identically the constraint equation (fl6l) . it follows that the fluid fields {p = p Q > 0,V,p,T} 
are necessarily solutions of the INSFE equations. This can be proven either: a) by direct 
substitution of / = /a/(x, t) in Eq.([T]) (Proposition A); b) by direct evaluation of the velocity 
moments of the same equation for G(x, t) = 1, v,m 2 /3 (Proposition B). In fact, thanks to 
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assumptions 1-6 the first two moment equation coincide respectively with the isochoricity 
and Navier-Stokes equations [Eqs. (jlj) and (jSJ)]. As a consequence the energy equation is 
also satisfied by the fluid fields V and p, namely there results identically in Q x / 



V 



^ V +-[V (p - 0) - f R ] - vV 2 Y 



0. (19) 



Therefore, the third moment equation delivers the Fourier equation [Eq.(l6"])]. The same 
proof (for Proposition B) is straightforward if / ^ /m(x, t). This is reached again imposing 
the same constraint equation ( fTBT) on first velocity-moment of the distribution function /. 
Finally, the proof of Proposition C is achieved by imposing on po(t) the constraint equation 
( I17p and invoking assumptions 1) and 5). It follows that the thermodynamic entropy can 
be identified with the Shannon entropy, i.e., letting for all t E I, Sr(t) = S (/m(x, t)) . As a 
final remark, we point out that by properly prescribing the initial and boundary conditions 
for the KDF, one can show that also the appropriate initial and boundary condition for the 
fluid fields can be satisfied (see related discussion in Ref. ^]). Furthermore, by construction, 
due to the constraint ([TBI) , /(x, t) = /(x, t)/p a is a velocity-space probability density, i.e., 
there results identically in Vt x J, J dv/(x,£) = 1. 

III. LAGRANGIAN FORMULATION OF IKT 

The results of the previous Section permit us to formulate in a straightforward way also 
the equivalent Lagrangian form of the inverse kinetic equation Q]. The Lagrangian formu- 
lation is achieved in two steps: a) by identifying a suitable dynamical system (here denoted 
as INSFE dynamical system), which determines uniquely the time-evolution of the kinetic 
probability density prescribed by IKT. Its flow defines a family of phase-space trajectories, 
here denoted as phase-space Lagrangian paths (LP's); b) by proper parametrization in terms 
of the LP's the KDF and the inverse kinetic equation, the explicit solution of the initial- value 
problem defined by the inverse kinetic equation ([1]) is determined. First, we notice that - 
in view of the THM - it is obvious that the LP's must be identified with the phase-space 
trajectories of a classical dynamical system x a — > x(£) = T t)to x D generated by the vector field 
X(x, t). Hence we shall assume that the initial- value problem ([3]), which is realized by the 



S 



equations 



|r(*) = v(t), 
f t v(t) = F(r(t),t;f), 



(20) 



r(t ) = r , 
v(t ) = v , 

defines a suitably smooth diffeomeorphism. Here, by construction: 

• denoting by x(i) = x( x o,^o,^) the solution of the initial- value problem (1201 . x„ = 
x(x(i),t, i ) is its inverse. Both are assumed to be suitably smooth functions of the 
relevant parameters; 

• both x(i) = x(x ,t ,t) and x Q = x(x(t), t, t Q ) identify admissible LP's of the dynam- 
ical; 

• r(t) is the Lagrangian trajectory which belongs to the fluid domain f2; 

• v(t) and F(r(t),t; /) are respectively the Lagrangian velocity and acceleration, both 
spanning the vector space R 3 . In particular, F(r(t), t; f), which is defined by Eqs. (I12j) - 
( fl5l) . and depends functionally on the kinetic probability density /(x, t), is the La- 
grangian acceleration which corresponds to an arbitrary kinetic probability density 

/(x,t); 

• /(x, t) is a particular solution of the inverse kinetic equation Q] which is subject to the 
assumptions imposed by the THM. 



It follows that the Jacobian [J(x(i), t) 



dx(t) 
<9x 



of the map x — > x(t), which is generated 



by Eq. (1201) for a generic distribution function /(x, t) of this type, reads 



J(x(t),t) 



exp I J dt'H(x(t'),t') 



(21) 



Jo 



where if(x.t) = f V Q - 4 [V-lfl Q + — u ■ V-IL Instead, in the case in 

v ' ' 2 pi pi * pf L =EJ ^ pi = ' 



which there results identically / = /m(x, t), the Jacobian reduces to 



J(x(t),t) 



dx{t) 



<9x 



<pW_ 

*4. P (*°) 



exp 



u 2 (t G ) w 2 (t) 



(22) 
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Second, it is immediate to prove that the kinetic equation in the Lagrangian representa- 
tion can be written in the form 



J(x(t),t)/( X (t),t) = /(x o ,t Q ) EE / C (X ) 



(23) 



where /(x(t),t) is the Eulerian representation of the KDF and / (x ) is a suitably smooth 
initial KDF. Eq. (l2~3~]) manifestly implies also the time evolution of /(x(t),£) in terms of the 
initial distribution function: 



From the mathematical standpoint main consequences of the theory are that: 1) the La- 
grangian formulation (of IKT) is uniquely specified by the proper definition of a suitable 
family of phase-space LP's; 2) Eq.( l24l) uniquely specifies the time-evolution of the Eulerian 
KDF, /(x(i),t), which is represented in terms of the initial distribution function / (x ) and 
the LP's defined by the INSFE dynamical system; 3) the time-evolution of the fluid fields 
{p = Po > 0, V,p > 0,T > 0,S T } is uniquely specified via the KDF /(x(t),t); 4) Eq.(jMJ) 
also provides the connection between Lagrangian and Eulerian viewpoints. In fact the Eule- 
rian KDF, /(x, t), is simply obtained from Eq.( |24l) by letting x = x(t) in the same equation. 
As a result, the Eulerian and Lagrangian formulations of IKT, and hence of the underlying 
moment (i.e., fluid) equations, are manifestly equivalent. From the physical viewpoint, it is 
worth mentioning the LP's here defined can be interpreted as phase-space trajectories of the 
particles of the fluid, to be considered as a set of "classical molecules", i.e., point particles 
with prescribed mass, which interact only via the action of a suitable mean-field force kind. 
The ensemble motion of these particles has been defined in such a way that it uniquely 
determines the time evolution both of the kinetic distributions functions and of the relevant 
fluid fields which characterize the thermofluid. 

IV. CONCLUDING REMARKS 

In this Note an inverse kinetic theory has been developed for the INSFE problem based 
on a restricted phase-space representation, i.e., in which the phase space of the kinetic 
description is identified with the direct product T = Q x V defined in such a way that Q 
coincides with the fluid. We have shown that equivalent IKT approaches can be formulated 



/(x(t),t) 



1 



f ( X (x(t),t,t )). 



(24) 



J(x(t),t) 
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both in the Eulerian and Lagrangian viewpoints. The result appears relevant for several 
reasons, in particular: 1) the Lagrangian dynamics here determined permits to advance in 
time self-consistently the fluid fields, i.e., in such a way that they satisfy identically the 
requires set of fluid equations. For isothermal fluids, this conclusion is consistent with the 
results indicated elsewhere by Tessarotto et a/.[5j; 2) the Lagrangian dynamics, defined in 
the configuration space of the fluid (fluid domain), takes into account the specific form of 
the phase-space distribution function which advances in time the same fluid fields; 3) the 
theory permits an exact description of the motion of the "classical molecules" which follow 
the Lagrangian dynamics. Finally, in our view the formulation here presented is promising 
in turbulence theory. In fact, the connection between Eulerian and Lagrangian phase-space 
descriptions based on a IKT approach, here established on rigorous grounds, represents a 
potentially useful new development. 
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